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IS. ‘ABSTRACT 

Seven  guidance  equation  techniques,  both  present  and  advanced,  are  dis¬ 
cussed.  Sufficient  detail  is  given  so  that  the  basic  philosophy  and  equation 
development  for  each  technique  may  be  understood.  The  techniques  described 
are:  Q,  Delta,  Explicit,  Linear  ^Tangent,  "Optimal,  "  Numerical  Integration, 
and  Parameter  Optimization.  Both  normal  and  fast  reaction  targeting  con¬ 
cepts  are  discussed.  Comparison  is  made  between  the  seven  techniques  with 
respect  to  computer  storage  and  speed  requirements,  targeting  requirements, 
and  mission  flexibility.  It  is  noted  that  Q  and  Delta  guidance  concepts  require 
the  least  computer  storage  but  also  offer  the  least  flexibility.  The  Explicit 
and  Linear  Tangent  techniques  require  exact  integration  of  approximate 
dynamics  equations.  Optimal  and  Numerical  Integration  guidance  both  require 
numerical  integration  of  more  exact  dynamics  equations.  Parameter  Optimi¬ 
zation  guidance  has  the  severest  computational  requirements,  but  offers  many 
significant  advantages.  It  is  hoped  that  this  review  will  motivate  new  pro¬ 
grams  now  in  the  planning  stages  to  consider  the  advanced  guidance  concepts. 
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FOREWORD 


This  report  was  prepared  in  support  of  the  Air  Force  Space  and 
Missile  Test  Center  (SAMTEC)  Inertially  Designated  Tracker  Study  and  is 
published  by  The  Aerospace  Corporation,  El  Segundo,  California,  under 
Air  Force  Contract  F-4701-72-C-0073. 

This  report,  which  documents  research  carried  out  from  January 
through  June  1972,  was  submitted  for  review  and  approval  on  8  November  1972 
to  Lt.  Col.  W.  J.  Lantier,  SAMSO(DYGM). 


Approved  by 


J.  R.  Allder,  Director 
Guidance  Dynamics  Subdivision 
Guidance  and  Control  Division 
Engineering  Science  Operations 


Publication  of  this  report  does  not  constitute  Air  Force  approval  of  the 
report's  findings  or  conclusions.  It  is  published  only  for  the  exchange  and 
stimulation  of  ideas. 


WARD  J.  LANTIER,  Lt  Col,  USAF 
Chief,  Missile  Guidance  and  Technology 
Deputy  for  Technology 
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I.  INTRODUCTION 


This  report  provides  a  discussion  of  seven  guidance  techniques.  Some 
of  these  techniques  are  considered  to  be  advanced  in  that  their  targeting  re¬ 
quirements  are  minimal.  It  is  believed  that  these  advanced  techniques  are 
potential  candidates  for  future  space  or  weapon  system  guidance  programs, 
and  that  guidance  software  for  advanced  vehicle  systems  should  not  be  restricted 
to  current  guidance  schemes.  The  advanced  guidance  techniques  included  in 
this  report  require  a  higher  speed,  greater  capability,  and  generally  a  more 
versatile  guidance  computer.  Crude  estimates  of  computer  speed  and  capacity 
are  provided  in  Table  1  for  each  of  the  techniques  discussed  in  the  report. 

The  seven  guidance  techniques  are: 

1.  Q 

2.  Delta  (or  A) 

3.  Explicit 

4.  Numerical  Integration 

5.  Linear  Tangent 

6.  "Optimal" 

7.  Parameter  Optimization 

A  summary  of  the  features  of  these  guidance  techniques  is  given  in  Table  1. 

Note  that  Q  guidance  is  included  for  completeness,  although  this  technique  is 
obsolete. 

It  has  been  noted  that  there  is  considerable  variation  in  the  computer 
requirements  for  each  of  these  techniques.  The  requirements  are  minimal 
for  Q  and  A  guidance;  in  fact,  these  techniques  were  developed  for  computers 
that  were  not  equipped  with  a  hardware  divide  instruction.  Their  computer 
technology  is  of  the  1956  era.  As  a  result,  these  techniques  were  extremely 
difficult  to  target  (i.  e. ,  to  obtain  all  guidance  constants). 


As  time  progressed,  computers  became  available  that  were  more 
versatile.  Then  the  "explicit"  guidance  techniques  became  practical  and 
were  devised  by  several  contractors.  Targeting  routines  were  still 
required  in  order  to  obtain  the  gravity-turn  pitch  program  for  flight  through 
the  atmosphere.  However,  a  complicated  least  squares  computation  pro¬ 
gram  to  obtain  polynomial  expansion  coefficients  was  not  required  as  it  was 
for  the  Q  and  A  mechanizations. 

The  computer  requirements  for  Numerical  Integration  guidance  are 
quite  large  (20,  000  words,  1  psec).  The  targeting  requirements  are  mini¬ 
mal  since  only  a  booster  pitch  program  is  required. 

The  Linear  Tangent  Guidance  Program  uses  an  approximate  integra¬ 
tion  of  the  rocket  equations  of  motion.  The  program  is  as  easy  to  target  as 
the  Numerical  Integration  technique,  since  a  multistage  capability  is 
included.  In  addition,  the  computer  requirements  are  not  much  greater  than 
those  for  explicit  guidance  programs. 

The  "Optimal"  Guidance  Program  is  based  on  a  calculus  of  variations 
solution  that  requires  machine  computation  (iteration  and  integration).  The 
computation  complexity  is  therefore  somewhere  between  the  two  previous 
techniques.  This  technique  also  requires  only  an  input  pitch  program  and  is 
therefore  easy  to  target. 

Parameter  Optimization  guidance  is  seif -targeting  and  also  has  the 
greatest  computation  requirements  on  speed  and  capacity.  It  is  equipped 
with  an  optimum  capability  for  guided  flight  through  the  atmosphere  and  has 
the  greatest  capability  for  fast  reaction  of  any  of  the  previous  techniques. 
With  the  ever  advancing  computer  technology,  it  is  possible  that  an  airborne 
digital  computer  will  one  day  be  designed  to  permit  usage  of  this  technique 
for  the  weapon  systems. 
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II.  PURPOSE  OF  GUIDANCE  EQUATIONS 


The  purpose  of  the  vehicle  guidance  system  is  to  utilize  certain 
measurements  of  the  vehicle's  state  and  to  alter  the  course  of  the  vehicle  in 
order  to  achieve  the  desired  end  conditions  in  the  presence  of  disturbances, 
which  include  but  are  not  necessarily  limited  to  the  following: 

•  Winds 

•  Non-standard  propulsion  parameters  (thrust,  specific 
impulse,  engine  misalignments,  centering,  etc.) 

•  Drag  and  li  „  uncertainties 

•  Weight  errors 

•  Control  system  dispersions 

The  basic  guidance  system  may  be  either  radio  or  inertial.  When  a  radio 
guidance  system  is  used,  the  measurements  consist  of  range,  aiimuth  and 
elevation  angle  (R,  A,  E).  In  some  cases,  rate  measurements  may  also  be 
obtained.  These  measurements  are  then  processed  in  a  Kalman  filter  to 
obtain  the  vehicle  navigation  data,  as  shown  in  Fig.  1.  All  computations 
shown  in  Fig.  1  are  performed  in  the  ground  guidance  computer.  The 
steering  commands  are  sent  to  the  vehicle  over  the  guidance  data  link, 
decoded,  and  fed  to  the  vehicle  control  system. 

When  an  inertial  system  is  employed,  the  sensed  measurements  are 
obtained  from  the  vehicle  accelerometers.  These  measurements  are  then 
processed  on  board  the  vehicle  in  the  Flight  Digital  Computer  as  shown  in 
Fig.  2.  Note  that  the  mechanization  shown  in  Fig.  2  assumes  that  the 
accelerometers  are  mounted  on  an  inertially  stabilized  platform.  The  navi¬ 
gation  loop  would  appear  different  when  a  strapdown  inertial  system  is 
employed. 

The  guidance  equations  blocks  for  either  the  radio  or  the  inertial 
methods  may  be  quite  similar.  The  role  of  guidance  is  (either  implicitly 
or  explicitly)  to  specify  the  corrective  action  such  that  the  actual  vehicle 
state  will  become  the  desired  vehicle  state  at  the  burnout  point  (or  at  orbital 
insertion). 
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blank 


Figure  i.  Guidance  Mechanization  Used  with  Radio  Guidance 


Figure  2.  Guidance  Mechanisation  Used  with  Inertial  Guidance  (IMU  Assumed) 


This  report  discusses  only  the  philosophy  of  the  guidance  equations. 

Some  insight  into  the  derivations  is  given.  In  general,  the  airborne  dig;. cal 
computer  used  with  inertial  guidance  systems  may  contain  many  functions  in 
addition  to  what  has  been  referred  to  as  explicit  guidance.  A  possible  list¬ 
ing  of  some  of  these  functions  and  the  computer  words  of  storage  and  timing 
estimates  are  given  in  Table  2. 

In  addition.  Appendix  A  is  included  to  show  the  detailed  equations 
associated  with  an  actual  guidance  computer  program.  The  bloc!'-  flow  diagrams 
apply  to  the  Radio  Guided  Gemini  launches  using  the  Titan  launch  vehicle. 

The  total  computer  size  (Burroughs  A1  at  ETR)  was  3200  words  including  an 
auxiliary  memory  unit.  Thus,  one  may  get  a  feeling  for  the  size,  complex¬ 
ity,  and  logic  loops  when  the  computer  permits  storage  in  excess  of  10,000 
words.  Some  filtering  is  included  in  the  time-to-go  and  steering  command 
computation  blocks;  however,  this  is  easy  to  recognize  and  may  be  ignored 
for  inertial  guidance  applications. 

Table  2.  Typical  Missile-Borne  Computer  Programs* 


Words 

Timing  (tps) 

Executive  Program 

500 

10,000 

Navigation 

500 

2,  500 

Explicit  Guidance 

1500 

2,  500 

Atmospheric  Steering 

200 

200 

Vehicle  Sequencing 

200 

100 

Coast  Equations 

200 

100 

Ground  C/0  and  Calibration 

2000 

Backup  Modes 

100 

100 

Digital  Flight  Control 

5000 

50,000 

Digital  Attitude  Control 

1000 

30,000 

11,200 

60,  000  Worst 
Path 

Estimate  by  J.  Shaul,  The  Aerospace  Corporation 
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III.  TARGETING 


An  important  figure  of  merit  of  any  set  of  guidance  equations  is  the  time 
required  to  target  the  equations.  At  the  present  time,  a  wide  disparity  exists 
for  the  various  techniques.  This  figure  includes  the  time  required  by  all  com¬ 
puters  involved  to  compute  all  input  constants  required  by  the  guidance  equa¬ 
tions  in  order  to  "fly"  the  mission.  It  also  includes  the  time  to  transmit, 
insert,  and  verify  these  constants  into  the  launch  guidance  computer.  Thus, 
it  is  clear  that  the  targeting  time  may  be  the  vehicle  flight  reaction  time.  Of 
course,  in  some  cases  it  is  possible  for  vehicle  or  guidance  hardware  flexi¬ 
bility  limitations  to  be  the  pacing  item  in  some  fast  reaction  launch  situations. 

The  normal  method  for  vehicle  targeting  (using  explicit  guidance  equa¬ 
tions)  is  to  use  a  large  high-speed  ground-based  computer  equipped  with  a 
multi-vehicle  simulation  (MVS)  program  to  generate  an  open-loop  trajectory. 

(A  simplified  block  diagram  of  a  vehicle  simulation  is  shown  in  Appendix  B. ) 

This  trajectory  is  tested  using  dispersion  runs  to  verify  acceptability. 

The  trajectory  is  then  used  to  generate  the  guidance  constants.  The  vehicle 
simulation  then  runs  a  closed-loop  trajectory  using  the  guidance  constants. 
Dispersion  runs  are  again  run  to  verify  that  the  trajectory  still  satisfies  the 
vehicle  constraints.  If  the  closed-loop  trajectory  fails  the  test,  a  recycle  of 
this  procedure  is  required.  A  block  diagram  of  this  procedure  is  shown  in 
Fig.  3.  Targeting  times  using  this  technique  vary  from  two  weeks  to  several 
months. 

The  targeting  procedure  may  be  speeded  up  by  using  prestored  pitch 
programs  and  autopilot  gains,  as  shown  in  Fig.  4.  These  are  generally 
functions  of  the  payload  weight,  altitude  of  injection,  and  desired  orbital 
inclination.  The  reaction  time  is  further  decreased  by  programming  the  actual 
vehicle  guidance  equations  into  the  vehicle  simulation  and  using  these  directly 
as  the  search  tool  to  obtain  the  closed-loop  reference  trajectory.  In  addition, 
a  print  denoting  satisfaction  of  the  constraints  should  be  furnished.  If  the 
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Figure  3.  General  Techniques  for  Parameter  Selection  (Targeting) 


PAYLOAD  WEIGHT 


«/> 


fr* 

O 


v 

Sh 

3 

tJ 

<u 

o 

o 

»H 

& 

bo 

a 

•H 

4-1 

V 

ao 

h 

rt 

H 

C 

o 

•  H 

•M 

o 

n) 

<u 

4-> 

m 

£ 


<u 

Ih 

3 

bo 


-11- 


reference  trajectory  is  not  satisfactory,  the  operator  must  modify  the  selection 
of  launch  azimuth  or  injection  altitude  manually  in  order  to  obtain  a  satisfactory 
nominal  trajectory.  In  the  event  that  the  fast  reaction  procedure  cannot  find  a 
satisfactory  reference  trajectory,  it  is  desirable  to  use  an  N- dimensional 
search  and  optimization  procedure.  Mission  reaction  times  of  several  hours 
may  be  obtained  when  this  technique  is  employed.  Even  faster  reaction  times, 
of  the  order  of  minutes,  are  achievable  when  dictated  by  program  requirements. 

The  passage  of  time  and  the  increase  in  comput  :  efficiency  has  seen 
guidance  equations  become  more  and  more  sophisticated  with  the  result  that 
targeting  requirements  have  become  minimized.  An  enormous  targeting 
burden  was  associated  with  the  Q  and  A  guidance  equation  formulations.  Each 
of  these  techniques  may  contain  30  to  well  over  several  hundred  constants  that 
must  be  determined  by  the  targeting  computer.  Even  the  gravity  computation 
contained  constants  that  required  specification  during  the  targeting  procedure. 

Either  the  general  purpose  or  the  fast  reaction  technique  may  be  used 
with  "explicit,"  "numerical  integration,"  "linear  tangent,"  and  "optimal" 
guidance  techniques.  The  fast  reaction  targeting  technique  works  particularly 
well  with  the  linear  tangent  guidance  equations  because  of  their  built-in  multi¬ 
stage  capability.  In  the  case  of  "Optimal"  and  "Explicit"  guidance  techniques, 
it  may  be  necessary  for  the  fast  reaction  targeting  routine  to  specify  "Aim 
Points"  for  each  of  the  exo-atmospheric  vehicle  stages. 
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IV.  Q- GUIDANCE 


Q-guidance  is  described  in  Refs.  1  and  2,  and  has  been  used 
successfully  in  Short-Range  Ballistic  Missile  applications. 

Existing  formulations  of  the  A  and  Q  guidance  equation  mechanizations 
involve  the  concept  of  the  required  velocity  vector.  The  required  velocity 
vector,  V  (x.t)  is  defined  as  the  velocity  that  is  required  at  the  present 
position,  x,  and  time,  t,  to  satisfy  the  mission  requirements  assuming  an 
impulse  burn  to  achieve  V  The  vector  x  is  measured  from  the  center  of 
the  earth  in  an  inertial  coordinate  system.  The  ballistic  missile  is  to  be 
steered  such  that  the  velocity  to  be  gained 


goes  to  zero  at  burnout  where  V  =  x  denotes  the  instantaneous  velocity  vec¬ 
tor.  With  this  steering  philosophy,  it  is  obvious  that  only  a  limited  number 
of  injection  constraints  may  be  satisfied.  For  example,  the  three  components 
of  the  desired  velocity  vector  at  burnout  can  be  achieved.  If  a  variable  coast 
period  to  the  target  is  allowed,  then  one  additional  constraint  may  be  satis¬ 
fied.  A  complete  description  of  the  constraints  that  may  be  satisfied  is  given 
in  Refs.  1  and  2. 

The  flexibility  offered  by  the  impulsive  burn  definition  of  the  required 
velocity  vector  is  completely  satisfactory  for  the  minimum  energy  type  of 
trajectory.  In  this  type  of  mission,  it  is  only  required  to  hit  a  target.  The 
dispersions  in  the  velocity  vector  when  the  target  is  hit  are  of  no  concern. 
This  approach  is  generally  not  sufficient  for  satellite  injection  purposes 
where  both  the  altitude  and  velocity  must  be  steered  to  the  desired  values. 

A  typical  altitude  dispersion  that  would  result  from  use  of  the  impulsive 
approach  is  several  miles. 
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The  guidance  equations  further  limit  the  specification  of  (x,t). 
This  vector  must  be  defined  such  that  any  freely  falling  particle  for  which 
V  =  VR  initially  will  continue  to  have  V  =  at  ail  subsequent  positions 
and  time. 

From  the  definition  of  Vp,  it  follows  in  general  that 


dt 


3 


£ 

R=  1 


a?" 


at 


(2) 


where 

1 

x  =  x 
2 

x  =  y 
3 

x  =  z 

and  x,  y,  z  denote  coordinates  in  the  earth  centered  inertial  computational 
system. 

For  the  particular  class  of  required  velocity  vectors  acceptable  for 
Q  guidance,  it  may  be  shown  that 


7? 


+ 

VR 


^R 

at 


The  rocket  equation  for  ballistic  missile  flight  is  written  as 

-*dV  — 

8  =  dT  "  aT 


(3) 


(4) 
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where  a^,  denotes  the  thrust  acceleration.  The  above  equations  may  be 
combined  to  give  the  Q-guidance  equation 


VE  '  ”aT 


8vr 

R=  1  ° 


(5) 


where 


V  ,  V2  =  V  ,  and  V3  =  Y 
gx  g  gy  g  gz 


Q-guidance  has  one  very  desirable  feature  in  that  it  does  not  require  a 

complete  navigation  determination.  That  is,  the  missile  does  not  need  to 

know  its  velocity  or  position  vectors  in  order  to  determine  V  .  This  quan- 

S 

tity  is  obtained  using  a  single  integration  of  the  Q-guidance  equation.  It  is 

necessary  to  specify  the  initial  conditions  V  (0)  and  the  Q  matrix 

g 


SV 


Qij  = 


R 


9xJ 


(i*j»  =  2,  3) 


(6) 


It  has  been  shown  that  these  quantities  may  be  chosen  as  constants  for  an  IRBM 

mission.  An  involved  targeting  search  routine  is  required  in  order  to  provide 

this  determination  and  still  keep  the  guidance  equation  injection  errors  within 

acceptable  bounds.  In  general,  the  Q..  matrix  coefficients  are  complicated 

functions  of  time  involving  many  empirical  constants.  These  constants  must 

be  determined  as  a  result  of  extensive  computer  programs. 

The  steering  procedure  used  with  some  Q  guidance  mechanizations 

requires  that  the  unit  thrust  vector,  £  ,  be  directed  along  V  : 

g 


V 

'z  _JL 

V 

g 


(7) 
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This  is  equivalent  to  issuing  a  commanded  turning  rate,  u>c,  to  torque  the 
rate  integrating  gyros  in  the  missile  control  system  as 

3  s  !*1  s  j  ( V  X  V  ) 
c  *  b  y  2  \  g  g/ 

g 

For  practical  reasons,  the  gain  on  the  vector  cross  product  may  be  replaced 

by  a  simple  time  varying  function  to  avoid  the  singularity  that  occurs  when 

V  =  0.  This  singularity  is  particularly  troublesome  when  the  thrust  mis- 
g 

alignments  are  significant. 

In  many  space  missions,  it  is  desirable  that  the  vehicle  contain  instan¬ 
taneous  position  and  velocity  information.  When  this  information  is  necessary 
for  other  purposes,  then  the  advantage  of  the  Q  mechanization  is  less  attractive. 
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V.  DELTA  GUIDANCE  EQUATIONS 


Delta  guidance  has  been  effectively  used  for  the  past  fifteen  years  and 
is  still  in  use  (Ref.  3).  It  is  a  viable  method  for  space  missions;  however, 
targeting  and  the  lack  of  one  flight  trajectory  control  variable  motivate 
against  its  use. 

In  this  mechanization,  the  vector  is  represented  by  expansions  of 
the  type 


VRx  =  VRxn  +  A<*  '  xn>  +  B(^  '  »n>  +  C(z  "  +  D(t  '  ‘n1  (9) 


plus  second  order  terms  as  required.  Similar  expansions  are  used  for  and 
Vrz’  Variations  in  this  specification  exist,  but  the  general  philosophy  is  the 
same. 


As  in  the  Q-guidance  problem,  the  question  of  constant  determination 
(A,  B,  C,  D, .  . . )  may  require  many  hours  on  a  large  capacity  computer.  The 
targeting  procedure  involves  the  determination  of  many  possible  burnout 
points  and  least  squares  fitting  of  this  data  to  the  selected  polynomial 
expansions. 

This  method  of  guidance  has  been  used  in  the  past  in  small  capacity 
airborne  computers  that  did  not  have  built-in  divide  and  square  root  instructions. 

The  A  mechanization  does  require  a  navigation  loop  to  determine  position 
and  velocity.  The  gravity  computation  is  also  an  expansion  with  one  formulation 
of  this  system.  That  is,  because  of  computer  limitations,  g  is  computed  using 
an  expansion  of  the  form 


gx  =  K(x  -  xQ)  [c( 


+  C1  x  + 


2  2  ~\ 
C2  z  +  C3  x  +  z  +  Cj  x  z 


(10) 
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with  similar  terms  for  g^.,  gz  rather  than  the  explicit  form 


8  =  •%  * 
r 


r  =  VxT/z 


Steering  with  the  A  guidance  mechanization  is  discussed  in  detail  in  Ref.  3. 
In  one  system,  the  unit  thrust  vector  is  directed  as 


v 

ar 


8  ar  J 


V  dt 
g 


where  the  *  on  V  denotes  that  the  pitch  (p)  and  yaw  (y)  components  of  V 

.  ®  O 

have  been  modified  by  adding  on  the  tangential  component,  as  follows: 


VgP  =  VgP  +  KlVgt 


V  +  K-V  . 

gy  2  gt 


where  and  are  determined  by  the  procedure  to  minimize  vehicle  propel¬ 
lant  losses  which  occur  during  the  closed- loop- g  ’dance  phase  of  flight.  Note 
that  and  V  still  vanish  at  final  sf.age  burnout,  since  this  event  is  com¬ 
manded  when  =  0.  The  integral  terms  in  the  steering  equation  have  been 
found  to  compensate  well  for  thrust  misalignments. 

The  A  guidance  equation  mechanization  has  been  tested  with  success 
for  flights  as  complex  as  the  synchronous  equatorial  mission.  The  procedure 
essentially  required  that  values  be  specified  to  hit  space  target  (position) 
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vectors.  Prediction  polynomials  were  the. a  used,  also  in  the  form  of 
expansions,  to  compute  velocity  components  at  the  end  of  coast  periods.  The 
prediction  polynomials  were  also  used  to  calculate  the  time  of  flight  during 
the  coast  periods.  Note  that  this  mechanisation  requires  that  verniers  be 
used  to  cut  off  each  stage  precisely  at  burnout.  This  follows  from  the  fact 
that  an  assumption  is  made  in  the  equations  that  the  targee  (position)  vector 
is  actually  hit.  That  is,  it  is  assumed  that  =  0  at  orbital  insertion. 
Obviously,  the  target  would  not  be  hit  if  a  large  velocity  error  were  made  in 
the  equations.  The  above  problem  is  probably  solvable,  but  the  machine  time 
and  complexity  of  equations  would  increase  several  times  in  the  event  that  a 
vernier  was  not  used  for  orbital  insertion,  especially  when  a  long  coast 
period  followed  main  engine  burnout.  This  problem  area  would  probaoly  be 
even  more  critical  and  difficult  to  cope  with  for  the  Q-guidance  mechanization. 
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VI.  EXPLICIT  GUIDANCE  EQUATIONS 


The  explicit  guidance  equation  mechanization  has  been  used  on  a 
number  of  space  missions.  This  technique  is  consistently  effective  and  will 
probably  be  used  for  several  more  years,  especially  with  all-inertial  systems. 
Details  on  explicit  guidance  are  given  in  Refs.  4  through  8.  The  derivation 
of  the  explicit  guidance  (and  steering)  procedure  is  generally  credited  to 
Duncan  MacPherson. 

The  principal  use  of  explicit  guidance  mechanizations  has  been  for 
guiding  a  single  rocket  stage.  When  additional  stages  are  used,  the  concept 
of  Aim  Points  has  been  used.  That  is,  the  equations  are  designed  to  cause 
each  stage  to  pass  through  a  pre -designated  aiming  point.  The  Aim  Points 
are  obtained  during  the  targeting  procedure. 

Explicit  guidance  equation  formulations  imply  that  closed  form  solu¬ 
tions  to  Newton's  equations  of  motion  are  employed  in  the  computations. 

At  Aerospace,  the  explicit  concept  implies  that  closed  form  solutions  are 
also  used  in  the  computation  of  commanded  thrust  attitude  or  thrust  attitude 
rate.  Thus,  tho  Aerospace  explicit  guidance  equation  mechanization  takes 
into  account  the  finite  thrusting  characteristics  of  engines  (as  opposed  to  the 
impulsive  burn  assumptions  of  the  Q  and  a  mechanizations).  This  pro¬ 
cedure  allows  additional  constraints  to  be  satisfied  at  burnout. 

Five  constraints  can  be  satisfied  with  the  Aerospace  explicit  guidance 
equation  mechanization.  These  include  not  only  the  three  components  of  the 
velocity  vector,  but  also  the  two  lateral  components  of  the  position  vector. 

The  tangential  component  of  vehicle  position  cannot  be  constrained  without 
excessive  loss  xr.  performance  because  the  vehicle  is  not  equipped  with 
engine  thrust  magnitude  control. 

The  explicit  solution  of  the  rocket  equations  is  computed  in  the  pitch 
plane  using  three  controls:  (1)  the  initial  (constant)  value  of  thrust  attitude. 
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(2)  the  thrust  attitude  rate,  and  (3)  the  final  (engine  cut-off)  time.  The 
lateral  position  and  velocity  are  also  controlled  by  using  corresponding 
attitude  degrees  of  freedom  in  the  yaw  plane.  A  summary  explanation  of 
the  explicit  guidance  technique  is  given  below  for  motion  in  the  pitch  plane. 

The  differential  equations  of  motion  for  the  thrusting  vehicle  in  the 
atmosphere  are  given  as: 


(r  6)  ’  =  raT  cosi|j+  (D,  L) 
r=  r02  -  ^=r+  aT  sini{j+  (D,  L) 


U') 


where 


r  =  distance  from  center  of  earth  to  missile 

0  =  angular  displacement  between  the  present 
position  to  a  future  position  measured  from 
the  center  of  the  earth 

a,p  =  thrust  acceleration 


(D,  h)  =  small  aerodynamic  drag  and  lift  terms 


=  thrust  attitude  measured  from  the  local 
horizontal 


Equations  (9-14)  may  be  written  in  terms  of  V,Y  where 

V  =  inertial  velocity  magnitude 

Y  =  elevation  angle  of  the  velocity  vector 

with  respect  to  the  local  horizontal 
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using  the  transformation 


r  =  V  sin  V 
r0  =  V  cos V 

The  above  equations  then  become 


V  =  a_  cos  ( 4*  -  Y  )  -  — y1  sin  y  +  (D,  L) 

1  r 

r  =  aT  sinvji  +  +  (D,  L)  (15) 

r 

These  equations  may  be  approximately  written  as 

V  =  aT  +  Small  Terms 

r  =  a,j4  +  Small  Terms  (16) 

It  has  been  shown  that  choosing  the  pitch  or  yaw  thrust  attitude  as  a 
linear  time  function  will  result  in  efficient  flight.  This  assumes  that  the 
thrust  vector  is  oriented  such  that  it  is  nearly  normal  (within  «  10  deg)  to 
the  local  gravity  vector.  Then 


^o  +  io  t 


(17) 


The  form  of  a,p  is 


‘T  - 


TO 


1  - 


aTOt 


(18) 
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Thus,  the  first  of  these  equations  may  be  readily  integrated  to  obtain 
an  approximate  time -to -go,  t  ,  until  the  stage  burns  out  and  the  desired 
value  of  the  final  velocity  is  achieved. 

When  the  approximate  value  of  t  is  obtained,  the  second  equation 
may  be  integrated,  neglecting  the  small  terms,  as 

‘t  *  *o  +  fe  aT(T>  [+0  +  *0  T]  dr 
*/  o 

(19) 

rf=ro+rotK+  f  -T)  aT(T)  +  A  -rldT 

Jo  L  o  O  J 

Equations  (9-19)  are  then  solved  for  ‘l'  and  such  that  the  mission  alti¬ 
tude  and  flight  path  angle  requirements  will  be  approximately  satisfied. 

The  effect  of  the  small  terms  is  then  evaluated  as  a  perturbation  solution 
and  a  small  adjustment  in  time-to-go,  attitude,  and  attitude  rate  is  made. 

The  use  of  the  Aerospace  explicit  guidance  equations  is  not  wasteful 
of  propellant  for  exo -atmospheric  flight.  This  is  evidenced  by  the  fact  that 
this  procedure  has  been  used  effectively  to  generate  nominal  trajectories 
at  Aerospace.  It  is  even  possible  to  develop  a  targeting  routine  by  flying 
the  first  stage  using  a  gravity  (zero  lift)  turn  through  the  atmosphere.  The 
explicit  guidance  technique  is  then  used  for  controlling  the  upper  stages. 

This  targeting  technique  requires  iteration  on  the  gravity  turn  and  also  on 
Aim  Points  if  several  upper  stages  are  required.  However,  propellant 
waste  is  not  a  serious  problem.  The  explicit  technique  allows  the  satisfac¬ 
tion  of  additional  constraints  as  required  for  satellite  injection  whefre  both 
the  lateral  position  and  velocity  must  be  constrained  at  burnout. 

A  second  principal  advantage  of  this  mechanization  arises  in  the  area 
of  targeting.  Only  a  minimum  number  of  trajectory  dependent  constants  is 
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required.  No  additional  constants  to  provide  optimum  trajectory  shaping 
need  be  determined. 

There  is  a  third  advantage  of  explicit  guidance  techniques.  This 
method  explicitly  solves  for  the  required  initial  attitude  and  attitude  rate 
required  to  accomplish  the  end  conditions  at  burnout.  The  Q  and  A  steering 
formulations  command  error  quantities  to  the  control  system  that  are  func¬ 
tions  of  V  .  Then  V  — -0  as  a  result  of  a  control  system  damping  procedure. 

3  S  - . 

That  is,  in  the  Q  and  A  mechanizations,  the  lateral  components  of  V  will  go 

S 

to  zero  as  the  solution  of  a  first,  second,  or  third  order  system  (whichever 
is  the  case),  exhibiting  the  transient  response  characteristics  of  these 
damped  systems.  Hence,  system  stability  is  of  primary  concern  with  sys¬ 
tems  of  this  nature.  This  type  of  solution  behavior  is  not  obtained  with  the 
highly  predictive  explicit  mechanizations.  In  this  mechanization,  the 
approach  to  the  desired  end  conditions  is  monotonic  without  overshoot,  under 
shoot,  or  any  significant  damped  oscillation.  This  follows  from  the  fact 
that  the  system  is  not  designed  to  null  small  departures  from  some  pre¬ 
specified  nominal  trajectory. 
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VII.  NUMERICAL  INTEGRATION  GUIDANCE 


This  guidance  technique  has  been  used  to  guide  large  rocket  vehicles 
using  radar  data  (R,  A,  E)  measurements.  As  a  result,  a  high-speed,  large 
(32  k)  digital  computer  is  employed  to  perform  all  the  computations. 

Attitude  profiles  are  devised  as  shown  in  Fig.  5  for  both  pitch  and 
yaw.  Note  that  there  are  three  controls  in  pitch,  9q,  Gq  and  0^  as  well  as  a 
similar  number  in  yaw  plus  a  final  engine  cut-off  control  T.  Hence,  seven 
parameters  at  orbital  insertion  can  be  steered  with  this  technique.  These 
include 

h  =  attitude 

V  =  flight  path  angle 

0  =  final  pitch  attitude 

P 

y  =  sideways  position 

y  =  sideways  velocity 

0^  =  final  yaw  attitude 

V  =  velocity  magnitude 

In  order  to  illustrate  the  operation  of  the  guidance  technique,  it  is  only 
necessary  to  consider  the  controls  and  constraints  in  the  pitch  plane. 

Hence,  we  assume  the  following  controls  and  corresponding  constraints: 
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It  is  now  convenient  to  consider  a  large  powered  flight  simulation,  with  ail 

powered  stages  contained  in  the  guidance  computer.  A  nominal  run  is  made 

•  • 

using  the  best  estimate  of  0q,  0 ,  0^,  and  T  in  the  computer.  Four  dispersed 

runs  are  then  made  in  which  each  of  the  four  control  parameters  is  varied 
one  at  a  time.  Note  that,  when  the  final  attitude  is  constrained,  there  is 
essentially  one  less  control  parameter  because  the  following  approximate 
relation  exists  between  0Q,  0^  and  T: 

0q  +  0q  (T  -  AT)  +  0f  AT  =  0f*  (20) 

Thus,  when  final  attitude  is  controlled,  the  following  three  equations  are 
obtained: 


9Vf 

90 

o 

AG  + 
o 

dA 

90£ 

AGf 

+ 

9 Vf 

9T 

9h{ 

-  A0  + 

9h^ 

A0f 

4. 

9hf 

90 

o 

o 

T 

9T 

9Yf 

A0  + 
o 

• 

A  A 

4- 

9^ 

90 

o 

90f 

4a  U£ 

T 

9T 

AT  =  AV 


D 


AT  =  Ah. 


D 


AT  =  AV, 


(21) 


(22) 


(23) 


where  the  partial  derivatives  are  obtained  by  numerical  differencing  of  the 
four  dispersed  runs  from  the  nominal.  The  deviations  on  the  right-hand 
side  of  the  equations  are  the  deviations  of  the  nominal  from  the  desired  burn 
out  conditions.  The  set  of  three  equations  is  solved  for  A0Q,  A0£,  and  AT. 
The  values  for  0q,  0q,  6^  and  T  are  then  obtained  using  the  nominal  values 
and  the  previous  attitude  constraint  relation. 

#In  practice,  an  exact  implementation  may  be  used  for  the  attitude  constraint 
relation. 
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Rate  steering  commands  in  pitch  are  then  sent  to  torque  the  vehicle 
control  system  rate  integrating  gyros  as  follows: 


CO 

c 


e-  ec 

At 


+  e. 


where 

0  =  actual  vehicle  pitch  altitude  estimated  by  the  Kalman  filter 

At  =  steering  command  computation  interval  («  0.  5  sec) 
k  =  gain  factor  chosen  to  provide  system  stability  (  k  »  1). 

A  similar  computation  is  made  for  yaw  steering  commands.  In  general,  gjc  will 
be  limited  to  some  value,  generally  less  than  2  or  3  deg/sec.  Thus,  when  the 
computed  value  of  a>c  exceeds  the  limited  value,  then  the  limited  value  is  used. 
When  T  -  AT  <  t  <  T,  it  follows  that  one  of  the  controls  is  lost.  Generally 
closed-loop  attitude  steering  may  be  terminated  at  this  time  and  the  constant 
values  for  0Q  and  6^  employed  on  an  open -loop  basis.  The  cut-off  computation 
to  obtain  updated  values  for  T  is  still  performed  approximately  every  0.  5  sec 
until  cut-off  countdown. 

Several  variations  of  this  general  guidance  technique  may  be  used. 

It  is  often  very  wasteful  of  propellant  to  steer  the  final  insertion  attitude 
as  a  control  variable.  In  addition,  it  is  often  unnecessary  for  mission  suc¬ 
cess  to  control  final  attitude.  In  many  missions,  the  payload  can  perform  its 
desired  function  regardless  of  the  insertion  attitude.  In  this  case  the 
numerical  integration  guidance  technique  may  be  employed  eliminating  0£  as 
a  control  (i.  e. ,  set  AT  =  0).  The  attitude  constraint  equation  will  not  be 
applicable  for  this  case.  The  three  controls  (in  the  pitch  plane)  are  then 


-30- 


T  — —V 


and  the  operation  proceeds  in  a  manner  similar  to  the  previous  description. 

The  advantage  of  numerical  integration  guidance  is  that  it  is  flexible  in 
adapting  to  vehicle  configuration  constraints.  The  disadvantage  is  the  require¬ 
ment  for  a  high-speed  large  capacity  computer. 
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VHI.  LINEAR  TANGENT  GUIDANCE 


Linear  tangent  guidance  was  devised  at  Aerospace  by  F.M.  Perkins 
and  is  described  in  Refs.  9  through  ii.  Although  the  technique  has  not 
been  employed  for  guidance  of  any  launch  vehicle,  it  has  been  used  in 
several  targeting  simulations.  One  problem  with  this  guidance  formula¬ 
tion  is  that  a  "flat-earth"  gravity  model  is  used  in  the  derivation.  Correc¬ 
tion  terms  for  this  approximation  are  subsequently  added  on.  A  princi¬ 
pal  advantage  of  these  equations  is  that  they  have  a  multi-stage  capability 
without  the  need  for  Aim  Points  as  is  generally  the  case  with  the  "explicit" 
guidance  equations. 

The  simple  derivation  for  the  basis  of  the  linear  tangent  formulation 
given  herein  differs  considerably  from  that  given  by  Perkins. 

It  may  be  assumed  that  it  is  desired  to  obtain  the  thrust  vector  steering 
profile  that  will  minimize  the  expended  thrust  velocity.  Then  a  suitable  cost 
function  is  obtained  as 


aT(r)dt.  (25) 

However,  the  purpose  of  guidance  is  to  steer  the  vehicle  such  that  the  posi¬ 
tion,  x*jy  velocity,  Xp,  and  possibly  vehicle  thrust  attitude,  Xp,  are 
constrained.  The  augmented  cost  function  may  then  be  formed  as 

J*  =  f  aT  dt  +  *  '  (*f  -  XQ)  +  ^  *  (xf  "  x’q)  +  (xf-xc) 

%/Q 

(26) 

where  A,  p,  rj  are  constant  Lagrange  multiplier  vectors,  and  the  subscript  f 
denotes  the  value  at  t  =  T.  The  flat  earth  rocket  dynamics  equations  are  then 
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(27) 


x  =  AT  5+  g 

where  g  is  assumed  to  be  constant.  This  equation  may  be  integrated 
to  obtain  position  and  velocity  as 


•  • 


x,  = 


XQ  +  gT  + 


/ 

•'ft 


aT  (r)  I  (T)dT 


(28) 


and 


x,  =  x  +  x  T  + 
f  o  o 


IgT2  +  /\T  (t)(T-t)|(t)  dr 


(29) 


The  result  of  substituting  into  the  augmented  cost  function  Eq.  (26)  is 


•  ■r 

Jo 


i  +  \  «  $(T)  +  U  .  4(t)  (T-r)  +n.  4(r)  5(t-t) 
+  fixed  terms 


dt  (30) 


where  6  (T-t)  denotes  the  Dirac  Delta  Function.  It  is  then  clear  that  the 
augmented  cost  function  is  minimized  by  choosing 

£(t)~X  +  (MT-t)+  H  6  (T-t)  (31) 

where  the  symbol  ~ denotes  parallelism.  Since  £  is  constrained  to  be  a  unit 
vector  in  the  direction  of  the  thrust  acceleration,  this  equation  must  be 
normalized. 

The  derived  form  of  the  vector  (j(t)  is  subject  to  some  interpretation. 
Note  that  the  guidance  constants  are  the  Lagrange  multipliers  X,  p,  and  *?. 
Further,  it  follows  that  thrust  velocity  is  minimized  by  using  constant  attitude 
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steering  when  it  is  desired  to  control  only  the  final  velocity  vector.  If  it  is 
also  desired  to  control  the  final  position,  thenM  ^  0  and  an  attitude  rate  is 
required  to  achieve  optimality.  It  also  follows  that,  in  the  event  that  the 
final  attitude  is  constrained,  the  vehicle  should  change  its  attitude  to  the 
desired  value  of  attitude  and  shut  off  engines  simultaneously.  Since  an 
attitude  rate  maneuver  such  as  this  would  generally  break  up  the  vehicle, 
several  seconds  (»10)  must  generally  be  allowed  for  this  maneuver.  Note 
that  this  was  done  in  the  case  with  numerical  integration  guidance,  where 
AT  »  10  seconds. 

It  is  easy  to  demonstrate  the  relation  of  the  above  derivation  to  the 
linear  tangent  guidance  formulation.  Consider  the  figure  below: 


The  figure  shows  the  gravity  vector  acting  down  along  the  unit  vector  d.  The 
thrust  attitude  vector  can  then  be  decomposed  along  the  local  horizontal  unit 
vector,  h,  and  d  as 


*h  +  ^h  *T_t)  ^  +  \d  + 


(T-t) 


N(t) 


(32) 
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where  N{t)  is  the  normalization  factor  required  to  cause  £  to  be  a  unit 
vector.  Then  it  follows  that 


d  +  ^ 

vn^ty 


(33) 


Thus,  it  is  shown  that  the  local  tangent  vector,  as  also  derived  in  Ref.  9,  is  a 
bilinear  function  of  time. 

A  principal  result  of  the  Perkins  linear  tangent  guidance  program  is 
that  the  vacuum  flight  rocket  equations  with  constant  gravity  are  exactly 
integrable.  Hence,  Perkins  includes  a  Newton-Raphson  technique  for  obtain¬ 
ing  the  Lagrange  multipliers  \  and  p.  The  values  for  the  partial  derivations 
are  also  included. 

The  advantages  of  the  linear  tangent  technique  are  many.  It  has  a 
multistage  capability,  provides  fast  reaction,  and  is  readily  targe^able..  The 
disadvantages  arise  from  the  approximate  computation  of  gravity.  However, 
Perkins  includes  a  technique  that  has  been  extensively  simulated  and  works 
well. 
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IX.  OPTIMAL  GUIDANCE 


A  guidance  technique  generally  referred  to  as  "Optimal  Guidance" 
is  described  in  Refs.  12  through  18.  Although  this  technique  has  not  been 
used  to  guide  tho  Tight  of  7,  rocket  vehicle,  the  program  has  been  simulated 
for  use  with  good  results . 

The  basis  for  this  guidance  technique  is  best  illustrated  by  the 
Pontryagin  Maximum  Principle.  A  brief  description  of  this  principle 
is  given  in  Appendix  C. 

Using  this  procedure,  the  augmented  cost  function  is  obtained  as 


=  <Mx(T),  T)  + 


(V  -  x)  +  g  (x)  +  aT  £ 


where  0  (x(T),  T)  denote  end  contraints  to  be  satisfied.  The  rocket  flight 
dynamics  constraints  are  used  as 


x  =  V 

V  =  g(x)  +  aTt 

with  xq,  Vq  given,  where  x  denotes  the  inertial  vector  from  the  center  of  the 
earth  to  the  rocket  vehicle.  As  before,  a^,  denotes  the  vehicle  thrust 
acceleration  magnitude  and  |  its  direction. 

The  Hamiltonian  is  then  obtained  as 
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It  is  now  necessary  to  choose  the  value  for  £  which  minimizes  the 
Hamiltonian.  It  is  clear  that  this  value  is  obtained  as 

»■  #  (36> 

When  this  value  for  |  is  substituted  into  the  Hamiltonian,  the  reduced 
Hamiltonian  H*  is  obtained  as 


H  =  1  +  X  1  V  +  P  •  g(x)  -aT(i 

Application  of  the  Pontryagin  equations  then  gives: 

9H^  _  _  p  •  9g(x) 

dx  9x 


(37) 


(38) 


and 


9H* 

9v 


-P  =\ 


Thus,  the  following  result  is  obtained: 


p  *9g(x) 
9x 


(39) 


(40) 


and  the  optimal  direction  of  the  thrust  vector  is  specified  as 


Note  that  in  general  there  are  six  parameters  to  be  specified,  e.  g. , 

a  and  ll 
ro  r  o 

except  that  these  reduce  to  five  when  the  constraint  that  §  be  a  unit  vector  is 
applied.  Note  that  this  is  exactly  the  number  of  injection  parameters  that  may 
be  specified  at  engine  burnout  (orbital  insertion)  when  final  vehicle  attitude  is 
not  controlled. 

In  the  case  where  g  is  approximated  as  a  constant  (as  in  linear  tangent 
guidance),  the  following  result  is  obtained: 

P=  0  (42) 


or 


£  ~  £0  +  i\>  t  (43) 

which  corresponds  exactly  to  the  result  obtained  in  the  previous  section. 
When  the  customary  round  earth  expression  for  gravity  is  used 


g  = 


-CM  -* 
r3 


(44) 


(45) 


where  pl,  i  =  1,2,3,  denote  the  three  components  of  the  vector"?,  xl  are  the 
components  of  the  vector  x,  and  6^  denotes  the  Knonecker  Delta. 


-39- 


The  above  program  has  been  developed  at  IBM  to  include  many  cases 
for  orbital  insertion  alternatives.  For  example,  velocity,  altitude,  and 
orbital  inclination  can  be  controlled  as  well  as  many  other  possibilities. 

The  advantages  of  this  approach  are  that  gravity  is  handled  directly  in 
an  optimal  fashion.  In  addition,  it  is  reported  that  the  program  size  and  speed 
requirements  permit  use  in  an  on-board  digital  computer. 

The  disadvantages  of  the  technique  are  that  it  is  difficult  to  include  the 
effects  of  atmospheric  (drag  and  lift)  forces  in  the  optimization.  In  addition, 
it  is  difficult  to  compensate  for  other  path  constraints  on  staging  angle  of 
attack,  heating  vehicle  rate  limiting,  etc.  Even  so,  the  technique  is  a  very 
powerful  one  and  may  find  wide  usage  in  future  rocket  flight  applications  as 
well  as  targeting  software  systems. 
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X.  PARAMETER  OPTIMIZATION  GUIDANCE 


This  technique  has  never  been  used  to  guide  a  launch  vehicle.  The 
computational  speed  requirements  exceed  those  of  the  other  techniques 
described  by  possibly  an  order  of  magnitude.  It  is  required  that  a  highly 
modeled  powered  flight  vehicle  simulation  be  run  one  or  more  times  each 
guidance  computation  cycle  (»lsec).  The  technique  described  herein 
has  been  used  effectively  as  a  general  vehicle  targeting  technique.  However, 
the  guidance  application  requires  a  program  that  will  run  in  real  time.  It  is 
possible  that  this  i«  achievable  at  the  present  time  for  radio  guidance  appli¬ 
cations  using  a  large  ground-based  guidance  computer.  However,  it  is  felt 
that  some  computational  speed  enhancement  will  be  required  before  the  tech¬ 
nique  can  be  used  with  inertial  systems. 

This  guidance  technique  employs  an  n  dimensional  search  and  optimiza¬ 
tion  procedure  to  determine  all  the  discrete  parameters  of  a  rate  steering 
profile  such  as  the  one  shown  in  Fig.  6.  A  three-stage  rocket  vehicle  is 
depicted  in  this  figure.  Note  that  ten  attitude  turning  rates  are  shown  (five  in 
Stage  I,  two  in  Stage  H,  <*nd  three  in  Stage  HI).  Experience  indicates  that 
this  is  an  adequate  number  of  parameters  to  implement  even  the  most  exotic 
mission  requirements.  In  addition,  the  final  engine  shutoff  time,  T^,  is 
shown  as  a  control.  The  problem  now  is  to  find  the  values  of  all  the  control 
rates  and  T^  such  that  all  vehicle  constraints  are  satisfied  and  a  cost  function 
(related  to  maximum  payload  or  minimum  propellant  wasted)  is  minimized. 

Vehicle  constraints  can  be  quite  general.  The  following  inequality 
constraints  are  typical: 

1.  Dynamic  pressure:  q  <  <lmax 

2.  Angle  of  attack:  a  <  a  at  q  and  a  <  a  at  Staging 

3.  Heating  indicator:  J  qV  <  Hq 

4.  Radar  and  telemetry  antenna  look  angle  constraints 

5.  Vehicle  turning  rates:  Iwl  <  w 
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In  addition,  the  mission  equality  constraints  may  be  listed: 

1.  Velocity  magnitude 

2.  Flight  path  angle 

3.  Velocity  azimuth 

4.  Attitude 

5.  Out -of -plane  position 

Constraints  on  the  final  attitude  may  be  added  if  required.  It  should  also  be 
noted  that  the  five  orbital  insertion  constraints  defined  above  could  be  written 
in  many  different  coordinate  systems,  such  as  ECI,  launch  pad  inertial,  or 
even  orbital  parameters  such  as  apogee,  perigee,  inclination  period,  etc. 

When  all  the  prespecified  vehicle  turning  rates  are  obtained  by  means 
of  the  general  n-dimensional  search  and  optimization  procedure  (NDSOP),  it 
is  only  necessary  to  command  the  vehicle  to  follow  the  specified  rates.  The 
time  to  shut  down  the  engines  is  obtained  as  T^. 

It  is  clear  that  the  key  to  a  successful  parameter  optimization  guidance 
program  depends  upon  the  computational  speed  associated  with  the  NDSOP. 
The  core  of  the  program  is  a  steamlined  vehicle  simulation  of  sufficient 
accuracy  for  guidance  purposes.  It  is  felt  that  the  simulation  can  be  simpli¬ 
fied  in  several  areas  and  still  be  accurate  enough  for  guidance.  For  example, 
it  is  probably  not  necessary  that  the  thrust  model  be  altitude  dependent.  In 
addition,  the  breaktimes  associated  with  the  discrete  rate  parameter  specifi¬ 
cation  may  be  fixed  rather  than  variable. 

The  operation  of  the  n-dimensional  search  and  optimization  procedure 
will  be  briefly  explained.  Consider  a  cost  function 

J  =  «>  (na)  (46) 

to  be  minimized,  where  ua  denotes  the  control  vector.  In  the  example,  there 
are  at  =  1 . 11  controls  denoting  .  .  .  ,  T^.  The  augmented  cost 
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function  is  formed  by  using  Lagrange  multipliers  to  adjoin  the  constraints 
as 

J*  (u)  =  <j>  (ua)  +  Ap  *  [gP  (ua)  - 

where  GP(ua)  and  Gp  denote  the  actual  and  desired  constraints.  Index  nota¬ 
tion  employing  the  summation  convention  for  twice  repeated  indices  will  be 
used.  Greek  indices  range  over  the  number  of  control  variables,  Arabic 
indices  range  over  the  number  of  equality  and  inequality  constraints.  Inequality 
constraints  that  are  violated  are  treated  as  active  equality  constraints.  Thus, 
the  number  of  constraints  may  change  from  iteration  to  iteration,  and  program 
control  logic  must  be  included  to  perform  this  bookkeeping  and  counting  function. 

It  is  assumed  that  an  initial  estimate  of  the  values  of  the  control  vari¬ 
ables,  u®  and  A?  is  available.  It  is  now  desired  to  findAua  and  AAp  such  that 
the  augmented  cost  function  will  be  minimized.  Then, 

J*  =  <!>  (u 01  +  Au")  +  (\p  +  AA  )  [aPfu®  +  Au“)  -  GP  1  (48) 

o  o  P  L  o  DJ 

This  expression  for  J*  may  be  expanded  to  obtain 

J*  =  J  +  Au“  (V  +  Ap  GP1  +  i  Aua  Au^  ft  +  \P  GP  1 

o  L  o  aJ  2  l °  W} 

+  AApAGP(u“)  +  GPA^Aua  (49) 

where  JQ  =  J*  (u®)  andAGP  =  gP(u^)  -  G^. 
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The  cost  function  may  be  minimized  with  respect  to  Au  by 
differentiating  to  obtain 


9J* 


a(Au“) 


=  Ma+  GaAXp  +  AuP  Map  =  0 


(50) 


where 


Ma  ”  ^a  + 

M<vp  =  4&P  +  ^oG«p 

The  cost  function  may  be  minimized  with  respect  to  the  Lagrange  multiplier 
by  differentiating  to  obtain 


9J* 

a(AXb) 


=  AGP  +  GjAua  =  0 


(51) 


fY  Q 

It  is  now  necessary  to  solve  these  two  equations  for  Au  and  AX1.  The 
solution  may  readily  be  shown  to  be 


=  [AGP  -  GPM^  Mpj  [oPM^Gq]'1 
AuP  =  -M^[Ma  +  GpA\P] 


(52) 


These  correction  values  are  then  added  to  Xq  and  u^,  respectively,  to  obtain 

o  o 

an  improved  trajectory.  It  is  felt  that  the  procedure  should  converge  within 
several  iterations  (<10)  to  obtain  an  optimum  trajectory  prior  to  launch.  It  is 
then  conjectured  that  a  single  computation  of  the  procedure  every  1  to  10  sec 
during  the  total  powered  flight  operation  would  be  adequate  for  guidance. 

(The  10-sec  computation  time  interval  is  probably  achievable  at  the  present 
time;  however,  the  1-sec  interval  is  not.  ) 
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The  parameter  optimization  approach  has  been  simulated  as  a  trajectory 
optimization  and  a  targeting  tool.  It  has  been  noted  that  a  very  accurate  com¬ 
putation  of  the  gradients,  4>a  and  G^  >  is  essential.  In  order  to  obtain  the 
desired  speed,  it  is  recommended  that  the  simulation  compute  the  gradients 
analytically  each  time  a  run  is  made  rather  than  to  numerically  difference 
the  runs.  It  has  also  been  shown  that  inaccurate  second  derivatives, 
and  G?ft,  are  quite  acceptable.  These  can  probably  be  precomputed  prior 
to  launch. 
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APPENDIX  C.  PONTRYAGIN'S  MAXIMUM  PRINCIPLE 
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Apply  boundary  conditions: 


